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I. INTRODUCTION 

In the early 1980 's it was discovered that the density functional exchange correlation 
potential has a derivative discontinuity when the particle number crosses and integer-i^. In 
this paper a relationship between the potentials of systems with different particle numbers is 
proven. It will be shown that the Hartree-plus-exchange-correlation potential for an integer 
A^-electron (A^ > 2) system differs by a constant from the corresponding potential of an 
(A^ — l)-electron system if the densities are determined from the same external potential. 
As a corollary it follows that the functional derivative of the independent particle kinetic 
energy functional of the A^ and (A^ — l)-electron systems also differ by a constant. 

II. PROOF 

In the adiabatic connection approacb^i"- of the constrained minimization formulation of 
density functional theory-^-^^ the Hamiltonian H'^ for a system of A^ electrons is given by 

i^^ = r^ + tK^ + ^;^,ext [pn] . (1) 

Atomic units, h = e = m = 1 are used throughout. T is the kinetic energy operator, 
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T'' = -IJ2'^1 (2) 



and jVee is a scaled electron-electron interaction. 



The the external potential 
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'^Af,ext 



[PN] = ^v2^ti[pN];ri), (4) 

is constructed to keep the charge density fixed at pat (r) ,the ground state charge density of 
the fully interacting system (7 = 1), for all values of the coupling constant 7. It has the 
forn>i^iii 

^Jxt(M ; r) = (1 - 7) VuxUpn]; r) 

+ vl ( [pn] ■,r)-v2{ [pn] ; r) +vl^t ( [pn] ; r) , (5) 



where vl^^{[pN]',r) = Vcxt (r) is the external potential at full coupling strength, 7 = 1, 
and fgxt([PAf] ; r) is non- interacting Kohn-Sham potential. The exchange plus Hartree 
potential"^'-'^fu2:([pAr]; I"), is independent of 7, while the correlation potential t'7([pAr];r) 
depends in the scaling parameter 7. 

The chemical potential 

f. = El- El^, (6) 

depends on the asymptotic decay of the charge densityi^i"— , and hence is independent of the 
coupling constant ■■^^^. In Eq. (|6]) E'lj_^ is the groundstate energy of the {N — l)-electron 
system with the same single-particle external potential f J^t ( [pn] ; r) as the A^-electron sys- 
tem: 
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Note that by construction of v]^^ {[Pn] ; i") , Eq- © Pn is independent of 7, but the ground- 
state density of the (A^ — l)-electron system p]^_i, is a function of 7. 
The correlation energy E] [p'lf_i\ is defined as^^ 



E2 K-i] = (^ 
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where 



\I^"7 ) is the Kohn-Sham (A^ — 1) independent particle groundstate wavefunction 



that yields the same density as the interacting (A^ — l)-electron system at coupling strength 
7. Since the correlation part of the kinetic energy is given by 
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and hence 
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^^ 



■7 

Pn-1 



7 

^iV-l 



^T 



Pn-1 



^Pl- 



v: 



N-l 
ee 



Pn-1 



(10) 



the derivative of E^ \_plj_i\ with respect to 7 can be expressed as 
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Upon adding and subtracting (c.c. stands for the complex conjugate of the previous term) 

d'j Pn 
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and utihzing the normahzation of the wavefunctions which imphes that 
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Eq. flTTl) becomes 
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is the sum of the exchange E^ [p]^_i] and mutual Coulomb interaction energy U [pj^_i\ of 
the (A^ — l)-electron system. The charge density p]^_i is a function of 7, therefore^^ 
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where 
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(17) 



is the sum of the exchange and Hartree potentials for the (A^ — l)-electron system. Using 



Eq. ([T6]), Eq. (Ej) and the fact that 
(Tn|) can be cast as 

E2 K-i] - T^ K-i] 
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\1' 7 ) yield the same density p]^_i, Eq. 
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(^Z ([P^] ; r) + -fVux {[pn] ; r) - -fVux {[Pn^i] ; r)) 



(18) 



From the definition of E'^ [p]^_i] , Eq. (jH]), the correlation energy can also be written as 



E] K-i] = T' [p^-i] - t° [pl_,] + 7Ke K-i] - ^Eux K-i] 
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where 
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It now follows that 
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where use was made of the relation 
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The last term in Eq. (122!) can be transformed as follows: 
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Here use was made Eqs. (171) and ( lT3l) . Combining Eqs. ( l22l) and f l2^ leads to 
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At the solution poinlr^ii^ the following equation has to be satisfied for the (A^ — 1) -electron 
system: 

-^^ — f^ + IVux ([p]^_i] ; r) + t;2' ([p]^-i] ; r) + vlJ\p^\ ; r) = /x]^_^, (26) 

"Pn~i l^J 

where /i]^_]^ is a constant, the chemical potential of the (A^ — l)-electrons moving in the 

external potential vJ^tdpN] ', r). Taking into account that the number of electrons is fixed at 

(A^ — 1) independent of 7, it follows that 



= 0, 



(27) 



where it is assumed that the order of the integration and derivative can be reversed. From 
Eqs. ([25]), (ESD and ([2] 



^^? [p'U 



Comparing Eqs. flT8|) and f l28l) shows that 

Q = I d'r?^^^ {vl^^ ([p^] ; r) - vZ^, ( K_,] ; r)) (29) 

where 

v:^A[p];r)=v2{[p];r) + ^v^A[p];r) (30) 

is the Hartree plus exchange-correlation potential for a system with density p and coupling 
strength 7. 

The charge density p]^_i (r) is a functional of the potential v2^t {[Pn] ] r) as can be seen 
from Eq. ([T])^^'^^. Therefore 
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Now 
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AT-l 



is the density response function of the (A^ — l)-particle system. From stability considerations 
Xat-i (r, r') is negative semi-definite and has one zero eigenvalue which corresponds to the 
invariance of the density when the potential is changed by a constantr^ii^. This implies that 



,3 ^P^-i(r) 



0= /c/^r^:i|^/(r) 



= - j d'r j d'r'f (r) xl^, (r, r') (^vU[PN];r') + |^t^Z([P^]; r')) (33) 

is only possible if / (r) = constant since VuxUpn]', r') + -§-''^2 HPn]', r') 7^ constant. This proves 
the main point of this paper: 



uxc 



xc ( [Pn] ; r) = vl^^ ( [pI_-^] ; r) + constant (34) 

when Pn and p^_i are determined by the same external potential w^xt {[pn] ; r) • 



III. PROOF OF COROLLARY 

Let the energy functional F^[pj^]^-^ be defined as 



Tlieni^'^ 



and 



n [Pn] = r° [pn] + lU [pm] + lE, [pn] + E: [pm] . (35) 



f^ + ^^L(M;r)=/i (36) 

OPN (r) 



'^^^^'''^-^U.elt(M;r) = ^],_, (37) 



constant (39) 



^p'n-i (r) 
From the last two equations we find that 

and therefore, using Eqs. ( l34ll and (J35ll . it follows that 

(5piv (r) 5p]^_i (r) 

IV. DISCUSSION AND SUMMARY 

The relationship in Eq. fl34|) is valid for N > 2. This follows from the step in Eq. fl28|) 
where the correlation energy of the A^ — 1 electron system is taken as non-zero. In the proof 
use is made of the A^ and A^ — 1 electron wave functions, hence the proof given here is valid 
for integer A^. The proof for non-integer values of the electrons and will be presented in 
another paper. 

In summary, it was shown that the Hartree-plus-exchange-correlation potential for an 
integer A^-electron (A^ > 2) system differs by a constant form the corresponding potential for 
an (A^ — l)-electron system if the densities are determined with the same external potential. 
As a corollary it was shown that the functional derivative of the independent particle kinetic 
energy functional of the A^ and (A^ — l)-electron systems also differ by a constant. 
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